New physics of grain boundaries in bcc metals from the atomic-level: molybdenum as 

a case study 



O 

o 
o 



o 



D. Ye§illetent, T.A. Ariast 
^Departments of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 
^Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, NY 14853 
^Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA 02139 

(February 1, 2008) 

We present a systematic trend study of the symmetric tilt grain boundaries about the (110) axis in molybde- 
num. Our results show that multiple structural phases, some incorporating vacancies, compete for the boundary 
ground state. We find that at low external stress vacancies prefer to bind to the boundaries in high concentra- 
tions, and moreover, that external stress drives structural phase transitions which correspond to switching the 
boundaries on and off as pipe-diffusion pathways for vacancies. Finally, we present physical arguments which 
indicate these phenomena are likely to occur in the other bcc transition metals as well. 
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I. INTRODUCTION 

Molybdenum, with its high melting point and rela- 
tively inert chemical nature is often considered for high- 
temperature structural applications, but its extreme brit- 
tleness limits its usefulness. Experimental studies [0J|] 
demonstrate that this brittleness is an intrinsic property 
of the material and largely unrelated to the presence of 
impurities. These studies suggest, moreover, that this 
brittleness arises from weak inter-granular cohesion along 
the grain boundaries In this work, we shed light on 
the microscopic physics of these boundaries by presenting 
a detailed, atomic-level trend study of the behavior and 
structure of their low energy phases and the transitions 
among these phases. 

This study reveals new physics in the interaction of the 
grain boundaries with vacancies. The traditional mech- 
anisms of interplay between vacancies and grain bound- 
aries include pipe diffusion of vacancies along the bound- 
ary [p^-p^, and absorption and emission of vacancies 
during continuous climb of primary or secondary dislo- 
cations [0-0 ■ We find, in addition, that boundary va- 
cancies prefer to collect together at high densities on 
the boundary plane. Our results also indicate that grain 
boundaries can either emit or absorb large concentrations 
of vacancies into the surrounding bulk while undergoing 
structural phase transitions under applied stress. While 
our focus in the present study is on the particular sys- 
tem of symmetric tilt boundaries around the (110) axis 
in molybdenum, which are known to dominate the re- 
crystallization texture of this material , we expect for 
reasons detailed below that these conclusions hold quite 
generally for symmetric tilt boundaries in bcc materials. 



II. PROCEDURE 

The heavy computational demands of full-blown ab 
initio electronic structure calculations Q and semi- 
empirical tight-binding models |pO|-p^ restrict their use 
to the study of relatively small systems and relatively few 



configurations. In order to understand complex processes 
such as fracture, dislocation migration and inter-granular 
cohesion, computationally more feasible empirical poten- 
tials must be used. Moriarty has developed such an em- 
pirical model based on a multi-ion interatomic potential 
developed from first principles generalized pseudopoten- 
tial theory [0. The resulting model generalized pseu- 
dopotential theory (MGPT) potential successfully pre- 
dicts the cohesive, structural, elastic, vibrational, ther- 
mal and melting properties of molybdenum as well 
as the ideal shear strength and self-interstitial and va- 
cancy formation energies ||2^. We use this potential 
throughout this work. 

Focusing on the (110) symmetric tilt boundaries, we 
consider 2=3(112) and 2=9(114), which are among the 
lowest in energy, and E=3(lll), S=9(221), E=ll(113) 
and S=ll(332) as examples of boundaries with higher 
energies. To study the physics of these boundaries, which 
reside in bulk material, we employ periodic boundary 
conditions as the most natural. To minimize boundary- 
image interactions, we always maintain at least seventeen 
layers of atoms between boundaries in our supercells. 

Determination of the ground state and the low energy 
excited state structures in principle requires the explo- 
ration of the phase space of all possible configurations, 
which is an impractical task without taking into account 
some basic physics. The primary consideration we use 
to restrict this phase space is that, due to the relatively 
strong directional bonding in molybdenum and similar 
bcc metals, the structure of the grain boundaries tends to 
preserve the internal topology of individual grains. Un- 
der this restriction, there remain then only three con- 
siderations for each boundary: (1) possible addition and 
removal of atoms to and from the faces of the grains at 
the boundary, (2) the displacement of the grains relative 
to one another, and (3) relaxation of the internal atomic 
coordinates. 

For the first consideration, the fact that the intersti- 
tial energy in molybdenum (> lOeV) is much larger than 
the vacancy energy (3eV) indicates that insertion 
of additional material at the boundary leads to unlikely 
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high interfacial energies. We therefore concentrate only 
on the removal of atoms at the boundary. Direct calcu- 
lations with the MGPT potential reveal that the most 
favorable sites for atom removal are in the vicinity of the 
boundary. It turns out that, because the atoms in the 
planes adjacent to the boundary plane (indicated by cir- 
cles in Figure |^a) pack closely together, these sites are 
the energetically most favorable for vacancy formation. 

This leads us to consider the following structural 
phases for the boundaries in our study: grains joined 
with the amount of material expected from the naive 
coincident site lattice (CSL) construction (^^full-material 
phase" ) , and boundaries where we remove atoms from 
the circled sites in Figure |l|a, which shows this construc- 
tion. Below, we find that binding energy per bound- 
ary vacancy is higher when vacancies collect together 
at high density on the boundaries. Therefore, we first 
concentrate on boundaries with high vacancy densities. 
Because removal of an entire plane of atoms near the 
boundary is topologically equivalent to the initial "full 
material" phase under appropriate relative displacement 
of the grains, we focus on the phase where we remove 
one half-plane of atoms from the layer adjacent to the 
boundary (^^ vacancy phase"). 



a) 



b) 



c) 
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FIG. 1. Structural phases of the E=9(114) grain boundary: 
(a) naive CSL boundary, (b) relaxed full-material boundary, 
(c) relaxed vacancy phase boundary. (To aid visualization, 
atoms from the two cubic sublattices are colored separately, 
light and dark.) 

To determine the ground state of the above two struc- 



tural phases, we next turn to the second consideration 
above, the relative shifts of the grains. The displacement- 
shift complete (DSC) cell, which we explore on a 16 x 16 
sampling grid (which reduces to 4 x 4 by symmetry), 
contains all possible unique relative planar shifts of the 
grains. For each such shift, we complete our search and 
address the third and final consideration by performing 
full relaxations of both the perpendicular expansion of 
the grains and their internal coordinates. This exten- 
sive survey requires force and energy calculations of ap- 
proximately 100,000 boundary configurations, and would 
be infeasible to carry out with electronic structure tech- 
niques. 

To confirm the effectiveness of this survey in identify- 
ing ground state structures, we repeat the above proce- 
dure with supercells in which we remove an entire plane 
of atoms from the boundary. For all six boundaries in 
our study, our procedure indeed identifies the appropri- 
ate shift to recover the initial, topologically equivalent 
ground state found for the full-material phase before the 
removal of the plane of atoms. 

III. RESULTS 
A. Low Energy Phases 

Our results reveal several general trends in the physics 
of the (110) tilt grain boundaries in molybdenum. As a 
specific example, consider the E=9(114) boundary, which 
Figure |l| shows. Panel (a) shows the naive CSL construc- 
tion of the full-material phase, whose ground state as 
identified through our procedure, appears in Figure 
We find that this phase lowers its energy through both a 
perpendicular expansion of the boundary and integran- 
ular shifts parallel to the boundary, both of which tend 
to increase the local volume for the closely packed atoms 
near the boundary plane, restoring them to a more bulk- 
like environment. 

Table | shows that the outward expansion is a general 
trend among all grain boundaries in our study and that 
shifts occur along the boundary in the direction perpen- 
dicular to the tilt axis (y— direction. Figure |l|) to allow 
for a more bulk- like local environment. We find no sig- 
nificant shifts along the tilt axis (z— direction. Figure |^) 
for this phase of any of the boundaries. Finally, the 
last column of the table gives the mechanical compliance 
{{1/k)) of each boundary, where k is determined from 
the quadratic response of the energy per unit area per 
boundary as the cell expands, which takes the form 

AU = AUo + ^k{x - Ax)^ + . . . (1) 

where x is the expansion of the cell, and Aa; and AUo 
are the relaxed perpendicular expansion and boundary 
energy given in the table, respectively. The table lists 
the difference between 1/k for the boundary and for the 
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same cell filled with bulk, as it is this difference which 
defines the response of the boundary, independent of the 
bulk content of the cell. 

Turning now to the vacancy phase. Figure ^ shows 
the results of our ground-state search for the E9(114) 
boundary as formed by removing the circled atoms from 
Figure We now find grain-shifting as well as local 
internal relaxations to create more bulk-like local atomic 
environments, and as Table |l^ summarizes, we again find 
this behavior for all boundaries in our study. The per- 
pendicular shift is always inward compared to the ground 
state of the full-material phase, so as to close the material 
void associated with the vacancies. Finally, we frequently 
find for the vacancy phase, even for the highly stable 
S3's, parallel shifts relative to the full-material phase 
which produce more natural bonding arrangements for 
the boundary to accommodate the vacancies. This ten- 
dency for accommodation is so strong as to induce the 
only shift along the tilt axis we observe in this study, for 
the vacancy phase of the Ell(332) boundary. 

Figure compares the ground-state boundary ener- 
gies of the various phases. In all cases, we find the full- 
material phase (black bars) to be lower in energy than 
the vacancy phase (hatched bars), although often not 
by far. As in the experimental case, the lowest ground 
state interfacial energy occurs for the naturally occurring 
twin E3 (112) boundary Moreover, apart from this 
twin, all remaining ground state boundary energies are 
fairly constant (within 25%), as also found in experiment 
(within 30% H). 




E9(114) £11(113) E3(112) E3(lll) £11 (332) E9(221) 

FIG. 2. Energetics of grain 

boundary phases: full-material phase (black bars), vacancy 
phase (hatched bars), bulk- vacancy phase (gray bars). 

Figure |^ presents another relevant comparison. To 
transform physically into the full-material phase, the va- 
cancy phase must first expel its vacancies into the sur- 
rounding bulk material. Thus, when considering transi- 
tions, the relevant comparison is between the vacancy 
phase and a third phase, which consists of the full- 
material phase plus the corresponding number of va- 



cancies in the surrounding bulk material (bulk-vacancy 
phase). (Table III provides the relevant information for 
the resulting bulk vacancies of this third phase.) Figure^ 
shows that, although creation of vacancies on the bound- 
ary always increases the boundary energy, the energy for 
creating the corresponding number of vacancies in the 
bulk is always higher. Our results therefore are consis- 
tent with the fact that the boundaries act as reservoirs 
for vacancies, as occurs during pipe diffusion. 

To verify, as mentioned above in Section II, that 
boundary vacancies indeed prefer to cluster together, we 
have also considered boundaries with nearly isolated va- 
cancies within our supercell approach. Table |^ presents 
energy, displacement and compliance results for bound- 
aries where the vacancy concentration (3.3%) is one fif- 
teenth that of the boundary-vacancy phase. Table ^ 
presents the analysis of this data. As the larger binding 
energies reflect, the concentrated phase is more stable. 
This added stability appears to arise from the structural 
relaxation through parallel shifting of the grains observed 
in Table |l|, made possible by the high density of vacan- 
cies. Finally, we note that in the extreme case of the Sll 
boundaries, vacancies do not even bind to the boundary 
at low densities. 



B. Phase Transitions 

Three trends which the preceding results exhibit can be 
expected from general, material-independent considera- 
tions: (1) that the boundaries present preferred binding 
sites for vacancies because they disrupt the bulk bond- 
order, (2) that boundary vacancies prefer to collect into 
the high-density vacancy phases because of the additional 
relaxation, which parallel shifting of the grains affords, 
and (3) that the binding of vacancies to the boundary 
reduces the intergranular spacing because this restores 
more bulk-like interatomic separations. These phenom- 
ena open the intriguing possibility that the application 
of tensile stress normal to grain boundaries in bcc metals 
generally drives transitions among these various struc- 
tural phases, thereby providing new forms of boundary- 
vacancy interaction. 

To explore this possibility, we consider the thermody- 
namic potential which is minimized under fixed external 
stress, the enthalpy. As a function of applied perpendic- 
ular stress (J, the enthalpy of a grain boundary structure 
relative to bulk is 



AH = AUo -aAx~ —A{l/k) + 0{a^) 



(2) 



Here AUo is the difference in the ground-state energy per 
unit area. Ax is the difference in preferred perpendicular 
intergranular separation, and A{l/k) is the difference in 
compliance. Each of these materials parameters appears 
for each relevant phase in Tables ^ ■ IV . 
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FIG. 3. Enthalpies of E3(112) boundary as a function ex- 
ternal stress for all three phases. 



into the bulk (cr^™**), thus making this transition more 
accessible. 

Finally, as we expect from the fact that boundary va- 
cancies are more stable in high concentrations, the stress 
required to induce formation of vacancies on the bound- 
aries {<J^c°"^) at low densities is generally greater than for 
the high-density vacancy phase. As a last consistency 
check on our analysis, we note that the breakdown stress 
for the bulk {a^J^"-^), as a characteristic of the perfect 
crystal and not the boundary, remains essentially un- 
changed between the two independent sets of calculations 
for low and high boundary vacancy concentrations. 



IV. CONCLUSIONS 



As an example, Figure ^ shows the behavior of the 
enthalpy, relative to bulk, of the naturally occurring 
I]3(112) boundary in its full-material, vacancy and bulk- 
vacancy phases. Three first-order phase transitions (en- 
thalpy crossings) are evident in the figure. At zero stress, 
as observed above, the full-material phase is the ground 
state of the boundary, and moreover, in the presence of 
vacancies, the vacancy phase has lower energy than the 
bulk-vacancy phase, indicating that vacancies prefer the 
boundary over the bulk. However, at an applied stress 
of about 18 GPa, the E3(112) boundary system under- 
goes a first-order phase transition in which the vacancy 
phase is no longer preferred, and the boundary ejects its 
vacancies into the surrounding bulk material. A second 
transition occurs near 25 GPa, at which point the bulk- 
vacancy phase becomes lower in enthalpy than the bulk 
material phase. This corresponds to the spontaneous for- 
mation of vacancies in bulk, indicating breakdown of the 
bulk material. The third transition occurs near 30 GPa. 
Were this transition accessible before the breakdown of 
the bulk material, it would correspond to spontaneous 
formation of boundary vacancies. 

Table VI presents the stresses for the above three tran- 
sitions for all boundaries in our study. In all cases, 
the emission stress is accessible before breakdown of the 
bulk material through spontaneous formation of vacan- 
cies. Moreover, except the for the outlying behavior of 
the E9(221) boundary, vacancies always first form spon- 
taneously in the bulk before they do so on the boundaries. 

To verify that these results are not artifacts of the high- 
density vacancy phase, we repeat the above enthalpy 
analysis for boundaries with low densities of vacancies 
using the data of Table IV. As Table VII summarizes, we 
again observe the same transitions. The Ell boundaries 
do not bind vacancies at low concentrations (Table 0), 
and therefore, the transition stresses for the emission of 
vacancies for these boundaries are not physically relevant. 
For the lower S boundaries (except the S9(221) bound- 
ary, which again exhibits an outlying behavior), diluting 
the vacancy concentration reduces the critical emission 
stress at which the boundaries emit boundary vacancies 



Several general conclusions reached in the discussion 
above are supported both by trends in our calculations 
for a specific set of boundaries in molybdenum and by 
quite general, material independent arguments. We ex- 
pect that the following conclusions likely hold gener- 
ally for the interactions among vacancies and tilt grain 
boundaries in bcc transition metals: (1) consistent with 
the traditional view of grain boundaries as diffusion path- 
ways, vacancies prefer the boundaries over the bulk at 
low stresses, (2) boundary vacancies prefer to collect into 
high-density vacancy phases, (3) application of sufficient 
tensile stress to a boundary induces a structural phase 
transition which drives the vacancies from the bound- 
aries into the bulk, thereby shutting off pipe diffusion 
along the boundary. The last of these conclusions in par- 
ticular may have important implications for crack growth 
through pipe-diffusion assisted void growth and void for- 
mation at grain boundaries. 

Finally, in terms of precise quantitative predictions of 
the critical stresses characterizing these phenomena, it is 
important to bear in mind that the particular interatomic 
potential which we have employed (MGPT), although 
one of the most reliable, is known to exaggerate energy 
scales for complex structures j2^]. We therefore would 
expect to find these same transitions, but most likely at 
lower stresses, when studied either experimentally or ah 
initio. 
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Boundary 


[mJ/m^] 


Ax 

[A] 


Ay 

[A] 


Az 

[A] 


A(l/fc) 
[mJ/mVPa^] 


E3(112) 


2130 


-0.1 


0.4 


0.0 


3030 


E3(lll) 


2450 


0.1 


0.3 


0.0 


1440 


E9(114) 


2360 


0.3 


0.1 


0.0 


7250 


E9(221) 


2380 


0.3 


0.9 


0.0 


2210 


Ell(113) 


2070 


0.3 


0.8 


0.0 


2310 


Ell(332) 


2610 


0.1 


0.4 


0.8 


1910 



TABLE II. Vacancy phase; energies (AWo), perpendicu- 
lar expansions (Ax), and shifts (Ay and Az) relative to the 
CSL construction, and compliances relative to bulk (A(l/fc)), 
where coordinates are as defined in Figure [ll. 



Boundary 


AUo 
[mJ/m^] 


Ax 

[A] 


A(l/fc) 
[mJ/mVPa^] 


E3(112) 


130 


-0.003 


1050 


E3(lll) 


100 


-0.002 


870 


E9(114) 


80 


0.003 


500 


E9(221) 


60 


-0.001 


470 


Ell(113) 


50 


-0.001 


130 


Ell(332) 


60 


-0.005 


180 



TABLE III. Isolated vacancy enthalpy information for each 
boundary orientation. Results are expressed for the num- 
ber of vacancies per unit area of the corresponding bound- 
ary-vacancy phase. Lattice expansion Ax and compliance 
A(l/fc) information is for longitudinal strain perpendicular 
to the boundary plane. 



Boundary 


AUo 
[mJ/m^] 


Ax 

[A] 


Ay 

[A] 


Az 

[A] 


A(l/fc) 
[mJ/mVPa^] 


E3(112) 


610 


0.1 


0.3 


0.0 


900 


E3(lll) 


2020 


0.5 


0.3 


0.0 


4770 


E9(114) 


1730 


0.5 


0.1 


0.0 


3410 


E9(221) 


2180 


0.4 


0.1 


0.0 


100 


Ell(113) 


1740 


0.5 


0.2 


0.0 


1460 


Ell(332) 


2160 


0.4 


0.5 


0.0 


580 



TABLE I. Full-material phase: energies (AUo), perpendic- 
ular expansions (Ax), and shifts (Ay and Az) relative to the 
CSL construction, and compliances relative to bulk (A(l/fc)), 
where coordinates are as defined in Figure 0. 



Boundary 


AWo 

[mJ/m^] 


Ax 

[A] 


Ay 

[A] 


Az 
[A] 


A(l/fc) 
[mJ/mVPa^] 


E3(112) 


720 


0.1 


0.3 


0.0 


950 


E3(lll) 


2076 


0.5 


0.3 


0.0 


2060 


E9(114) 


2360 


0.4 


0.1 


0.0 


3520 


E9(221) 


2200 


0.3 


0.1 


0.0 


260 


Ell(113) 


1900 


0.5 


0.2 


0.0 


1900 


Ell(332) 


2250 


0.4 


0.5 


0.0 


990 



TABLE IV. Bulk- vacancy phase: energies [AUo), perpen- 
dicular expansions (Ax), and shifts (Ay and Az) relative 
to the CSL construction, and compliances relative to bulk 
(A(l/fc)), where coordinates are as defined in Figure m. 
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Boundary 


Binding energy 


Binding energy 




per vacancy at 


per vacancy at 




high density 


low density 




[eV] 


[eV] 


S3(112) 


0.7 


0.6 


E3(lll) 


2.1 


1.1 


E9(114) 


1.4 


0.7 


E9(221) 


2.3 


2.0 


Ell(113) 


1.4 


-7.1 


Ell(332) 


1.4 


-1.7 



TABLE V. Boundary- vacancy binding energies at high and 
low densities. (The Ell boundaries do not bind vacancies at 
low concentrations.) 



Boundary 






^tear 




[GPa] 


[GPa] 


[GPa] 


E3(112) 


18 


25 


30 


E3(lll) 


17 


24 


not observed 


E9(114) 


21 


23 


25 


E9(221) 


38 


24 


10 


Ell(113) 


21 


28 


42 


Ell(332) 


31 


33 


35 



TABLE VI. Critical stresses for the the phase transitions 
discussed in the text: emission of vacancies from the boundary 
into the bulk (aj;'"''), breakdown of the bulk through sponta- 
neous formation of vacancies (erf formation of vacancies, 
and thus tearing, at the boundary (a*°'""). 



Boundary 




^break 


^tear 




[GPa] 


[GPa] 


[GPa] 


E3(112) 


16 


25 


32 


E3(lll) 


12 


24 


not observed 


E9(114) 


9 


23 


58 


E9(221) 


44 


24 


8 


Ell(113) 


90 


28 


60 


Ell(332) 


125 


33 


52 



TABLE VIL Critical stresses for the same transitions as in 



Table VL but at 3.3% vacancy concentration. 
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